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Abstract 

The rare decay — )• K*^i~^(!~ is a flavour changing neutral current decay with 
a high sensitivity to physics beyond the Standard Model. Nearly all theoretical 
predictions and all experimental measurements so far have assumed a K*^ P-wave 
that decays into the K'^tt^ final state. In this paper the addition of an S-wave 
within the K'^tt" system of — t- K*^i~^£~ and the subsequent impact of this on 
the angular distribution of the final state particles is explored. The inclusion of the 
S-wave causes a distinction between the values of the angular observables obtained 
from counting experiments and those obtained from fits to the angular distribution. 
The effect of a non-zero S-wave on an angular analysis of B^ — t- K*^£'^i~ is assessed 
as a function of dataset size and the relative size of the S-wave amplitude. An S-wave 
contribution, equivalent to what is measured in B^ — )• J/ijjK*^ at BaBar, leads to a 
significant bias on the angular observables for datasets of above 200 signal decays. 
Any future experimental analysis of the K^tt^ I'^i^ final state will have to take the 
S-wave contribution into account. 
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1 Introduction 



The description of flavour physics in the Standard Model (SM) has so far accurately 
matched the observations in the data from the B factories, the Tevatron and the LHC very 
well. However, there are several fundamental questions which do not have an explanation 
within the SM such as the mass hierarchy of the quarks and why there are three generations. 
To avoid creating large flavour changing neutral currents, any physics beyond the SM 
that contains new degrees of freedom that couple to the flavour sector is required to be at 
an energy scale of multiple TeV or to have small couplings between the generations, i.e. 
couplings that closely mimic those of the SM. The measurement of the inclusive 6—)- 57 
width [1| is one of the strongest constraints on new physics from the flavour sector; for the 
exclusive decays, K*^i'^i~ is of major importance. 

The analysis of B^ — )■ K*^i~^i~ is based on the evaluating the angular distribution 
of the daughter particles |2|. How to extract the maximal amount of information from 
the decay while keeping uncertainties from QCD minimal has recently attracted much 
interest [3]-[8]. The results from the experimental analyses of B^ ^ K*^i~^i~ [9 12 have 



focused on the forward backward asymmetry of the dimuon system {A^b) and the fraction 
of longitudinal polarisation of the K*^ (Fl) as a function of the dimuon invariant mass. 

With the acquisition of large data sets of B^ — > K*^i'^i~ decays, scrutiny is required of 
assumptions that have been made in current experiments. Nearly all theoretical papers to 
date use the narrow width assumption for the K~^7t~ system meaning that the natural 
width of the K*'^(892) is ignored. This means there is no interference with other i^+vr" 
resonances. Existing B^ — )■ K*'^i^i~ analyses consider B^ — t- K*^i'^i~ signal with K^tt' 
candidates in a narrow mass window around the K*^{892). However, in this region there 
is evidence of a broad S-wave below the i^'*'^(892) and higher mass states which decay 
strongly to K+tt', such as the S-wave K^^{U30) and the D-wave Kf{U30) [l3|. The 
best understanding of the low mass S-wave contribution comes from the analysis of K~^7i~ 



scattering at the LASS experiment 14 



The interference of an S-wave in a predominantly P-wave system has previously been 
used to disambiguate otherwise equivalent solutions for the value of the CP- violating phase 
in 5° ||15j and B^ 16 oscillations. In the determination of (ps in the -B^ — >■ J/ip(f) decay it 
was also shown that it is required to take the S-wave contribution into account [TT] and 



this has subsequently been done for the experimental measurements 18 20 



In this paper, we detail how a generic K~^7r~ S-wave contribution to B^ ^ K*^£~^£~ 
can be included in the angular analysis. Firstly, we develop the formalism set out in [21j 
to explicitly include a spin-0 S-wave and a spin-1 P-wave state in the )■ K~^7r~i~^i~ 
angular distribution. Here K*^ is used for any neutral kaon state which decays to K~^7i~. 
The impact of an S-wave contribution on the determination of the theoretical observables 
is evaluated in two ways: in the first we look for the minimum sample size in which an 
S-wave contribution (such as measured in [Ts]) significantly biases the angular observables; 
secondly we determine, for a given sample size, the minimum S-wave contribution needed 
to bias the angular observables. We then demonstrate how the S-wave contribution can 
be correctly taken into account and evaluate the effect of this on the statistical precision 
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that can be obtained on the angular observables with a given number of signal events. 



2 The K*^i+i- angular distribution 

The differential angular distribution for — > K*^i~^i~ is expressed as a function of the 
five kinematic variables (cos 6';, cosOk, 0, p"^ and g^). The angle Ok is defined as the angle 
between the and the momentum vector in the rest frame of the K*^. The angle 6i is 
similarly defined between the fi~^ in the rest frame of the dilepton pair and the momentum 
vector of the B^. The angle is defined as the signed angle between the planes, in the rest 
frame of the B^, formed by the dilepton pair and the K^tt^ pair respectively!^ The mass 
squared of the K'^tt^ system is denoted p"^ and the mass squared of the dilepton pair g^. 
The angular distribution is given as a function of cos 6i , cos 6k and as 

d,Mp^dcoXdco.9,d^ =i {''^ + "'^ + + 2/1) co= 29, + 2,, s.n^ 9, co= 2^ 

+ 2V2I4 sin 26i cos + 2^2/5 sin 61 cos + 2/6 cos 61 (2.1) 
+ 2^2/7 sin 9i sin + 2a/2/8 sin 29i sin + 2^2/9 sin^ 9i sin 20 

Ignoring scalar and tensor contributions, the complete set of angular terms are 
= I AlP + \Aor\' + S^^f^ Mlo^«o) + 44lAI' ' 
n = I + \^L±\' + (L -> R)) (^1 - + (Ax^ii± + All All) , 

^2 = \f^f + \-^L±\' + lAiiT + \-AR±n , 

h = l^f {\A^L±\' - I Alii P + |A±P - IAhH , 

h = (3fJ(AoA||) + iL^ R)) , (2.2) 

h = V2/3i (g?(AoA±) -{L^ R)) , 

h = 2/3/ (3ft(A|| A±) -{L^ R)) , 
J7 = V2l3i (53(AoA||) - ^ ^)) , 

Is = (53(AoA±) + ^ R)) , 
Jg = /3f (S(A|| A±) + iL-^ R)) , 



^This is the same sign convention for cos 61 and cos 9k as used by the BaBar, Belle, CDF and LHCb 
experiments [9- 12 and the same convention as used in LHCb (22^. 



where AH(o,\\,±,t) are the hehcity amphtudes and (3f = ^/l — ^mf /q^ j2j. In this paper 
the lepton mass is assumed to be insignificant, such that the angular terms with rn^/q^ 
dependence can be neglected and /?/ = 1 such that Ji and I2 can be related by = —II 
and /| = 

For a K^TT~ state which is a combination of different spin states, the amplitudes for a 
given handedness {H = L, R) can be expressed as a sum over the resonances (J) 

^ (2-3) 

j>i 

where Yj^{6k, 0) are the spherical harmonics, M is the matrix element and Pj{p^) is the 
propagator of the spin state which encompasses the dependence. A detailed description 
of the spin-dependent matrix elements and normalisation factors can be found in Ref. |21) . 

3 Angular distribution of K^TT~t^£~ for a com- 

bined S- and P-wave 

For K^TT^ masses below 1200 MeVj^ the contribution to the amplitudes from the D-wave 
K*^(1430) is so small that it can be ignored and only the J = 0, 1 terms in the sums of 
Eq. |2 .31 will be considered. The S-wave contribution to these amplitudes only enters in ^0 
giving 

V47r V 47r 



Ah,\\ = y ^^1,^^,11 cos 9k, (3.1) 
Ah,± = \ -^Ai^H,± sin 9k, 

V OTT 

where the spherical harmonics are expanded, leaving the propagator and the matrix element 
as part of the spin-dependent amplitudes 

Ao,Hfi oc Mo,Hfi{q^) Po(p'), 

Ai,H,^ cc M.^nAq^) PiiP 
A,,H,\\ cc M^^H,\\{q^) Pi(p2), 

where the first index denotes the spin and the normalisation from the three-body phase 
space factor is omitted. The propagator for the P-wave is described by a relativistic 



(3-2) 



^Natural units are assumed throughout this paper 
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Breit-Wigner distribution with the amphtude given by 



where mj^*Q is the resonant mass and 



r.,(/) ^ r^.,. ^-J ^ (3.4) 

the running width. Here t is the momentum in the rest frame of the K~^Tr~ system 



and to is t evaluated at the K^n pole mass. B is the Blatt-Weisskopf damping factor 23 



with a radius Rp. The amplitude can be defined in terms of a phase (6) through the 
substitution 

2 2 

cot (5= (3.5) 

to give the polar form of the relativistic Breit-Wigner propagator 

Piip') = (3.6) 
cot — t 

The LASS parametrisation of the S-wave |l4] can be used to describe a generic K~^'k~ 
S-wave. In this parametrisation, the S-wave propagator is defined as 

Poip') = + e'''H—^)p' (3.7) 

t cot Ob — i cot Or — I 

where the first term is an empirical term from inelastic scattering and the second term is 
the resonant contribution with a phase factor to retain unitarity. The first phase factor is 
defined as 

cot5B = ^ + ^rt, (3.8) 

where r and a are free parameters and t is defined previously, while the second phase 
factor describes the i^'o''(1430) through 

2 2 

Here, ms is the S-wave pole mass and Fs is the running width using the pole mass of the 
ii'o°(1430). The overall strong phase shift between the results from the LASS scattering 
experiment and measured values for J/ip K^tt" has been found to be consistent 

with TT [Is]. The parameters for the spectrum used in this paper are given in Table |ll 
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Table 1: Parameters of the K^ir resonances used to generate toy data sets. The K* 
masses and widths are taken from Ref. 13 and the Kf^ Blatt-Weisskopf radius and the 



LASS parameters are taken from Ref. 24 



State. 


mass 


r 


R 


r 


a 






(MeV) 


(MeV) 


(GeV)-i 


(GeV)-i 


(GeV)-i 






894.94 ±0.22 


48.7 ±0.8 


3.0 








-"■0 


1425 ± 50 


270 ± 80 


1.0 


1.94 


1.73 


TT 



The angular terms modified by the inclusion of the S-wave are /i,2,4,5,7,8 and the 
complete set of angular terms expressed in terms of the spin-dependent amplitudes is 

1 3 \/3 

I{ = —\AoLo\^ + — \AiLo\^ cos^ 9k + 2—\AoLo\\AiLo\ cos 6^ ^ cos 9k + {L ^ R) 
47r 47r 47r 

It = ~ (IAliiP ± |AiL±P + {L-> R)) sin' 9k 
1 



-'2 



— r 

-'15 



h = — 



h = ~ - IAliiP ± ^ R)) Sin' 9k 

2^(^0L0^ii||)cos5^nsin^;^ 

+ ^^^{A,loAIl\\) sin 9k cos 9k + {L ^ R) 
^^^{AoloAIl±) cos5o _l sin 9k 

+ ^^'^(AiloAIlx) sin 9k cos 9k - (L ^ R) 
h = 2|- {^{A,L\\Al^^) - (L ^ R)) sin' 9k 
^yi^(^oLo^^Lll) cos5^_ll sin^;^ 

± ^\I]^'^{AiloA\l\\) sin 9k cos 9k - {L ^ R) 



2 3 1 



1 

72 



(3.10) 



l7= — 



1 

V2 



5 



^y|^(^0L0^iL±) cos(5^_^ sin Ok 
/g = 1^ {%A^L\\A\^^) + (L -> i?)) sin^ 



The interference term of /i shows how this parametrisation encompasses the strong phase 
difference between the S and P-wave state. The left handed part of the interference term 
for /i can be written as 

2|AoLoPiLo|cos(5o^o«2|Mo,L,o||^'o(p')||Mi,i,o||Pi(p')|cos(5oV (S-H) 

where 



+ 5po - 5 



Mil 



-6 



Pi- 



(3.12) 



where Smjlq is the phase of the longitudinal matrix element and 6pj is the phase of the 
propagator. The phases in the interference terms for 14,5,7.8 can be similarly defined. For 
real matrix elements, i.e. nearly true in the Standard Model, the phases are equal for both 
handed interference terms 6^ = 6^. The phase difference between the S-wave and the 
P-wave propagators can be expressed as a single strong phase, 5s- 

The spectrum for the 5°— )■ K~^n~i~^i~ angular distribution can be calculated by 
summing over the S- and P- waves and integrating out the cos6'i, cos^i^- and dependence. 
This is illustrated in Fig. [l] where the matrix elements from Refs |4||6] at a value of 
6 GeV^ are used. Here the S-wave amplitude is assumed to be equivalent to the longitudinal 
P-wave amplitude. The S-wave fraction in the 800 < p < 1000 MeV window around around 
the P-wave is calculated to be 16% when using this approximation. As will be seen later 
there are no interference terms left in the angular distribution after the integral over 
cos 9k- 



4 The effect on K+n-i+i- observables 

So far the forward-backward asymmetry (Aps), the fraction of the K*^ longitudinal 
polarisation (Fl) and two combinations of the transverse amplitudes {A^ and Aim) have 
been measured. As such, these are the observables that will be concentrated on here, 
is defined in terms of the amplitudes as 

, .2, 3 mmAlL,_)-miR\\Aln±) .... 

AFB[q j = 7; 14 12 I M 12 I M — 12 — ^'^■^ 

2 \Aio\^ + \Ai\\\^ + \Aij_\'' 

for a pure P-wave state where the generic combination of amplitudes AjiA*j^ is defined as 

AjiA*j, = AjiLA*j,L + AjiRA*j,j^. (4.2) 
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p2 / GeV^ 



Figure 1: An illustration of the spectrum for the P-wave (dashed) and the S-wave 
(dash-dotted). The total distribution from both states is the solid line. The values were 
calculated at = 6 GeV^ by integrating out the angular distribution of — )• K^n^l'^(.~ 
using equal matrix elements for each state. The S-wave fraction here is 16% between 
800 < p < 1000 MeV 



where i G {0, ||,_L,t} and J = 0,1. The factorisation of the amplitudes into matrix 
elements and the propagators removes the dependence from the theoretical observables. 
In a similar way, Fl, A\ and Ai^ are defined as 

i^L(g^) 
AW) 







lAo 


2 


lAol 


2 + 


All 


2 + |A±|2 


IA± 


2 _ 


IA| 


2 


IA± 


|2 + 


lAii 


2 



:i-Ft 



\A 



10 1 







2 _ 


All 


2 


lAo 
(A|p 


2 + 


Aiil 


2 + 


A±|2 


lAii 


2 + 




2 ■ 



(4.3) 



These theoretical observables are normalised to the sum of the spin-1 amplitudes. In 
terms of the angular distribution, A-pB can also be expressed as the difference between the 
number of 'forward-going' /i+ and the number of 'backward-going' yU"^ in the rest frame of 
the 5°, 



-1. 



d cos 9i 



dr 



dr 



dg2dcos6'/ dg2 



(4.4) 



which explains the name of the observable. In Ref |22j, this expression was used to 
determine the zero-crossing point of A-pB- The inclusion of the S-wave in the complete 
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angular distribution means that Ap^ can no longer be determined by experimentally 
counting the number of events with forward-going and backward-going leptons, as Eqs. |4.1 



and AA are no longer equivalent. However, as the S-wave has no forward-backward 
asymmetry, no bias occurs in the determination of the zero-crossing point by ignoring the 
S-wave. The total normalisation for the angular distribution changes to the sum of S- and 
P-wave amplitudes, 

r" = = l^ior + l^mr + + \Aoo\'. (4.5) 

such that there is a factor of 

^''^"'"^ = ( l^oP+|A.||P + |A.J^+|^o„p j f^'^' 

between the pure P-wave and the admixture of the S and the P-wave. This is the fraction 
of the yield coming from the P-wave at a given value of and q^. Similarly, the S-wave 
fraction is defined as 

Vl^iol + l^iiir + I^i-lI + l^ool / 

and the interference between the S-wave and the P-wave as 

. / 2 f \ AoLo\\AiLo\ cos 5l + {L ^ R) \ , . 



2 

Substituting the above observables into the angular terms gives 

P 1 3 3 

= —J^s + -rJ^pFh cos^ 9k + — cos 9k, 

§r = ll:FAl-F,)il-cos^9K), 

P f 1 3 3 

-| = - — -T^s + -^J^P (1 - -^l) cos^ 9k + -4,s cos 9k cos 9k 

I \47r Ati Ait 

II 1 3 



J-p(1-Fl)(1-cos2^,^), (4.9) 



JpA| (1 -cos^ 9k) 



r" 4 Svr 

r" 2 Svr 

^ = 2~TpAfb {I -cos^ 9k) 
1 Stt 3 

= —J^pAi^ (l - COS^ 9k) . 
1 Stt ^ ^ 



For the purpose of this paper, a simplification of the angular distribution can be achieved 



by folding the distribution in (p such that = — tt for < 25 . The ^4,5, 7,8 angular 
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terms which are dependent on cos0 or sin0 are cancelled, leaving /i, 2,3,6) in the angular 
distribution: 

+2/6 COS + 272/9 sin 6*; sin 20 ) . (4-10) 
Combining Equation 4.10 with 4.9 gives the differential decay distribution, 



d^r 



9 



r" dg2dp^dcos0/<dcos^/d(/)' IGvr 



4.10| with |49] gives 

2 4 \ 
- J's + 3 A cos J (1 - cos^ 

2Flcos2 Ok {I — cos ^/) 
^(l-FL)(l-cos2^i,)(l + cos2^0 

+ ^(1 - - cos^ - cos^ Oi) COS ! 

4 

+ -v4fb(1 - cos^ Ok) cos9i 
+ Ai^{l - cos^6'i^)(l - cos^ 61) sin! 



3" " ' 3 
+ Jp 

+ 



The angular 
in Eq. |4lT 



distribution as a function of cos 6'; and cos 6k is given 



(4.11) 
be integrating over 



1 d^r 

r" dg^dp^d cos Oxd cos 6*, 



9 / /2 4 \ 

Piel ( g-^S+g^SCOS^^j (l-COS^^O 

2Flcos^ -cos^ 9i) 



+ ^(1-Fl)(1-cos2 M(1 + cos^ 
4 

+ -Afb(1 - cos^ 6'x) cos 61 



(4.12) 



^0 



and further integration from Equation 4.11 yields the angular distribution for each of the 
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angles, 



3 

+ -{1- Fi^){l + C0S^9i)+AFBC0Sei , 
o 

1 d^r _ ^Qg^ 
r" dg2dp2dcos6'ii- 2 ^ ^ ''^ 



^Fl cos^ + ^(1 - Fl)(1 - cos^ Ok) 



1 d^r _ 1 / 3 
r dg2dp2d0' ~ 4^ ^ 



i^L + ^ (1 - i^L)^T COS 20' + Ai^ sin 20' 



(4.13) 



The angular distribution can be integrated over using the weighted integral 

O(g^) = ^ , (4-14) 

J dp2dg2dp 

for the value of the observables integrated over a given region in "p^ . This leads to the 
integrated observables Fp, Fs and As which are solely dependant on . By definition, the 
fraction of the S-wave and the P-wave sum to one, Fs + Fp = 1. The complete angular 
distribution without any dependence is given by 



dT 



r' dg^dcos^i^rdcos^, 



2Flcos^ ^if(l -cos^ Qi) 
^(1-Fl)(1-cos2^k)(1 + cos2^^ 



+ (1 - i^s) 
+ 



^(1 - Fl)>1t(1 - cos^ - cos^ Qi) cos 20' 



4 

+ -Afb(1 - cos^ Qk) cos 6'/ 
+ - cos^ Qk){^ - cos^ 6'/) sin 20' 



where the 



normalisation of the angular distribution is given by 

r' = ^ 

d?2 



(4.15) 



(4.16) 
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P-wave 
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COS 01, 



P-wave 
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CO 



3 



re 




(c) 



Figure 2: One-dimensional projections of (a) cosOi, (b) cosOk, (c) 0' for the angular 
distribution of 5° — )■ with (blue-dashed) and without (red-solid) an S-wave 

component of 7%. The dilution effect of the S-wave on the asymmetry in cos 6'/ and the 
asymmetric effect in cos 6^ can be clearly seen. 



The 'dilution' effect of the S-wave can clearly be seen from the factor of (1--Fs) that appears 



in front of the observables in Eq. 4.15 



The effect of an S-wave on the angular distribution as a function of cos 9k, cos 9i 
and (j)' as illustrated in Figure [2} Here it is possible to see that the asymmetry in cos 9i, 
given by ApB, has decreased and that there is an asymmetry in cos6'^' introduced by the 
interference term. 



5 Effect of an S-wave on the angular analysis 

In an angular analysis of — > K~^iT~i'^i~ , the S-wave can be considered to be a systematic 
effect that could bias the results of the angular observables. The implications of this 
systematic effect are tested by generating toy Monte Carlo experiments and fitting the 
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Table 2: Parameters used to generate toy datasets. Ap-Q 
~ GeV^ 



Ref. 



22 



in the 1 < g < 6 



bin. 



The Fs value is taken from Ref. 



and Ai^ are taken from 
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Obs. 


^FB 


Fl 




Aim 




Value 


-0.18 


0.66 


0.294 


0.07 


0.07 



angular distribution to them. The results of the fit to the observables are evaluated for 
multiple toy datasets. 

The effect of the S-wave is evaluated for two different cases. Firstly, the effect of 
S-wave interference is examined as a function of the size of the dataset used. The aim 
of this is to give an idea of the current situation and the possible implications on future 
measurements of — )■ K~^7i'~i~^i~ . Datasets of sizes between 50 and 1000 events are 
tested. For comparison, the latest results from LHCb |22j have between 20 and 200 signal 
events in the 6 different bins considered. Secondly, the effect of different levels of S-wave 
contribution is examined. At present, the only information about the S-wave fraction is 



the measurement of Fs of approximately 7% in the decay — )■ J/tpK^ir^ from 15 for 
the range 800 < p < 1000 MeV. As the value may be different in i?° — )■ K^Ti~i^i~ , we 
consider values of Fs in this region ranging from 1% to 40%. The fraction of the S-wave, 
Fs, is expected to have some q"^ dependence because of the dependence of the transverse 
P-wave amplitudes. 

The parameters used to generate the toy datasets are summarised in Tables [T] and |2j 
The values of the angular observables used to generate toy Monte Carlo simulations are 
taken from the LHCb angular analysis of 5°—)- K*^fj,~^fj,~ in the 1 < g2 < 6GeV^ bin 
Within errors, these measurements are compatible with the Standard Model prediction for 
B^ —7- K*^[^l~ and the central value of the measurement is used. The nominal magnitude 
and phase difference of the S-wave contribution are taken from the angular analysis of B^ 



J/ip K+TT- 15 



The toy datasets are generated as a function of the cos6';, cos^i^, and p'^ using the 
angular distribution given in Eq. |4.11 For each set of input parameters 1000 toy datasets 



were generated. For each of these toy datasets, an unbinned log likelihood fit is performed 
that returns the best fit value of the observables and an estimate of their error. The 
expected experimental resolution is obtained by plotting the best fit values of an observable 
for the ensemble of toy simulations as illustrated for ApB in Fig. |3] (left) The pull value for 
an observable (O) is defined as 

Pb = ^^^^ (5.1) 
where ctq is the estimated error on the fit to the observable O*. This distribution is seen 



in Fig. |3] (right). The mean and the width are extracted from a Gaussian fit. For a well 
performing fit without bias, the pull distribution should have zero mean and unit width. 
A negative pull value implies that the result is underestimated and a positive pull value 
implies overestimation of the true observable. 
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100- 
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Apg value 



2 4 

pull 



Figure 3: Distribution of (left) the Ape results and (right) pull values for fits to 1000 
toy simulations each containing 1000 events. The S-wave is ignored in these fits. The 
resolution obtained is 0.026 ± 0.001. Since the S-wave is ignored there is a non-zero pull 
mean at (0.84 ± 0.03)o" . The width of the pull distribution is consistent with unity at 
(0.99±0.02)a. 
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Figure 4: Resolution (left) and pull mean (right) of 1000 toy datasets analysed as a pure 
P-wave state as a function of dataset size. It can be seen that the bias on the observable 
increases dramatically as the sample size increases. This is because the statistical error 
decreases increasing the sensitivity to the S-wave contribution. The bias of Afb is positive 
because Aps in negative in the bin chosen. 



The impact of ignoring the S-wave in an angular analysis of 

K*^£+£- 



5.1 



Firstly, the effect of an S-wave was tested as a function of dataset size in order to find a 
minimum dataset at which the bias from the S-wave in the angular observables becomes 
significant. Datasets were generated for sample sizes ranging from 50 and 1000 events and 
analysed assuming a pure P-wave state. The results are shown in Fig. |4| 
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o Apg ■ A^p 



Im 




Figure 5: Mean of the pull distribution of 1000 toy datasets analysed as a pure P-wave 
state as a function of S-wave contribution. The bias can be seen to increase with the size 
of the S-wave contribution in a linear fashion. 



From Eq. 4.12, it can be seen that has a factor of (1--Fl) in front of it. The large 
value of Fl used in generated the datasets is in turn causing to have a much worse 
resolution than Ap^, Fi^ and Aim- There is significant bias (non-zero mean) of the pull 
distribution for all observables when the S-wave is ignored for datasets of more than 200 
events. This corresponds to a change of 0.2cr in Fl for a dataset of 200 events. The 
behaviour can be understood in terms of the (1 — Fs) factor in Eq 4.12 It gives an offset 
to the fitted value of the observables which are proportional to the value of Fg. 

Secondly, the angular fit was performed on toy datasets with an increasing S-wave 
contribution. Datasets of 500 events were generated with a varying S-wave contribution 
in the narrow mass window of (800 < p < 1000 MeV) from no S-wave up to a Fs 
value of 0.6. The resolution, the mean and width of the pull distribution for each of the 
four observables {Aps, -^l, ^im) were calculated and the results are shown in Fig. [s] 
Significant bias is seen in the angular observable for an S-wave magnitude of greater than 
5%. The linear increase in the bias is another consequence of the (1-Fs) factor. 



5.2 Measuring the S-wave in B° — )- K+tt £'^£ 

Obtaining unbiased values for the angular observables beyond the limits shown requires 
a measurement of the S-wave contribution rather than ignoring it. With the formalism 
developed in Sect. 4, three options are explored for measuring this. The first option is to 
ignore the dependence and simply fit for p^-averaged values of Fs and As- The second 
option is to fit the line-shape simultaneously with the angular distribution. This can 
be done in a small p window between 800 and 1000 MeV or in the region from the lower 
kinematic threshold to 1200 MeV. In all cases the datasets used to perform the studies 
are identical to those used in Sect. 5.1. The difference is in how the fit is performed. In 
each case, the dataset and the S-wave sizes refer to the number of events in the smaller 
window. 
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Figure 6: Resolutions for three different methods to incorporate the S-wave relative to the 
resolution obtained when the S-wave is ignored. It can be seen that the best resolution is 
obtained when using the largest window. The original resolution is recovered to within 
10%. 



The angular distribution without dependence is given in Eq. 4.15 for each set of 
samples, we look at the resolution, the mean and the width of the pull distribution of the 
angular observables. 

The change in the resolution obtained on the angular observables for the three methods 
of including the S-wave in the angular distribution is demonstrated by plotting the ratio 
with respect to the resolution obtained when a single P-wave state is assumed. 

The resolutions and the mean of the pull distributions for the three different fit methods 
(ignoring the dependence, fitting a narrow window and fitting a wide window) 
relative to the resolution and mean obtained using the assumption of a pure P-wave 
state. The ratio between the fit methods including the S-wave in angular distribution and 
assuming a P-wave state as a function of dataset size are shown in Fig [6j The pull mean 
for all four fit methods is shown in Fig [7j 
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Figure 7: Pull mean for the three different methods to incorporate the S-wave and when 
the S-wave is ignored. There is a slight bias when the S-wave is included for datasets of 
less than 200 events but this bias is removed from all the observables when the S-wave is 
included in the fit for datasets of over 500 events. 



For all observables, it can be seen that the resolution degrades when the S-wave is 
included and the dependence is ignored. The resolution degrades by a smaller amount 
when the dependence is included in a small bin and the original resolution is recovered 
to within 10% when using the large range. There are two effects contributing to the 
improvement of the resolution. There are more P-wave events in the larger range and 
the wider mass window allows for the S-wave to be constrained by using the information 
from above and below the P-wave resonance. This results in the best resolution when the 
S-wave is included in the angular distribution. 

For all the observables, the pull mean approaches zero for datasets of greater than 300 
events implying that the bias present in all the observables when a pure P-wave state is 
removed when an S-wave is included in the angular distribution. This means that the 
inclusion of the S-wave component will be mandatory for all future experimental analyses. 
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6 Conclusion 



In summary, the inclusion of a resonant K^t{~ S-wave in the angular analysis of 
K*^i^i~ has been formalised and the complete angular distribution for both an S and P 
wave state described. We find that the inclusion of an S-wave state has an overall dilution 
effect on the theoretical observables. The impact of an S-wave on an angular analysis 
is evaluated using toy Monte Carlo datasets. We find that the S-wave contribution can 
only be ignored for datasets of less than 200 events. The bias on the angular observables 
incurred by assuming a pure P-wave K^ti~ state can be removed by including the S-wave 
in the angular distribution. The degradation in resolution on the angular observables from 
fitting a more complicated angular distribution can be minimised by performing the fit in 
a wide region around the K*^{%92) resonance. 
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